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Optimal link adaption to the scattering function of wide sense stationary uncorrected scattering (WSSUS) mobile 
communication channels is still an unsolved problem despite its importance for next-generation system design. In 
multicarrier transmission such link adaption is performed by pulse shaping which in turn is equivalent to precoding 
with respect to the second order channel statistics. In the present framework a translation of the precoder optimization 
problem into an optimization problem over trace class operators is used [1], [2]. This problem which is also well- 
known in the context of quantum information theory is unsolved in general due to its non-convex nature. However 
in very low dimensions the problem formulation reveals an additional analytic structure which admits the solution 
to the optimal precoder and multiplexing scheme. Hence, in this contribution the analytic solution of the problem 
for the 2x2 doubly-dispersive WSSUS channel is presented. 



1 Introduction 

It is well known, that channel information at the trans- 
mitter increases link capacity. However, since future 
mobile communication is expected to operate in fast 
varying channels, it is not realistic to assume perfect 
channel knowledge at the transmitter. On the other 
hand statistical information can be used which does 
not change in the rapid manner as the channel itself. 
In multicarrier communications this can be employed 
for the design of transmitter and receiver pulse shapes. 
However, the problem of optimal signaling in this 
context is still an unsolved problem. 

In this paper the most basic case of precoder and 
equalizer optimization with respect to the statistics of 
a doubly-dispersive channel is considered. Hence, the 
focus is on a model in which two random complex sym- 
bols (iid and zero mean distributed) are transmitted par- 
allel over a 2 X 2 random channel. Before transmission 
a special kind of linear precoding is performed which 
is motivated from Weyl-Heisenberg signaling scheme. 
The channel itself is a randomly weighted superposition 
of four possible channel operations, which are 

1) do not change anything 

2) permutation in the time domain 

3) permutation in the frequency domain 

4) permutation in the time and frequency domain 

Each action is distributed independently from the others 
and scaled with a certain amount of power. What is the 
optimal precoding and multiplexing scheme? It will be 
shown that this is the formulation of the WSSUS pulse 
shaping problem for Weyl-Heisenberg signaling in the 
lowest possible dimension. 

The paper is organized as follows. In the first part 
Weyl-Heisenberg (or Gabor) signaling in £2(1^) and 
is introduced. Then the main optimization func- 



tional according to [1], [2] is established. It will open 
up the relation to the pure state channel fidelity opti- 
mization which is an ongoing research topic in quantum 
information theory. In the main part the problem is 
solved for C^. 



2 Signal Model 

2.1 Weyl-Heisenberg Signaling on £2(1^) 

In this article a transmit baseband signal s{t) is con- 
sidered which is a superposition of time-frequency 
translates of a single prototype function 7(t) with 
II7II2 = 1. The translations are according to a subset of 
a lattice AF'^ in the time-frequency plane, i.e. 



S{t) = ^X„{SAnl){t) 



(1) 



where in this constellation F = Z and A denotes the 
2x2 generator matrix. The indices n = {ni, 712) range 
over the doubly-countable set T C F^, referring to the 
data burst to be transmitted. Let 



(2) 



denote the time-frequency shift operator (or phase 
space displacement operator). It is well-known that the 
operators S'^ establish an unitary representation of the 
polarized Heisenberg group and up to phase factors they 
are equal to the Weyl operators (see [3]). The complex 
data symbols to transmit are x„ where ni is the time 
instant and 712 is the subcarrier index. The transmit 
signal is passed through a linear time-variant channel 
denoted by the operator Ti, and further distorted by an 
additive white Gaussian noise process n{t). Hence, the 



received signal is then 



r{t) - ins)it) + n{t) = j i:{^i){S^s){t)dti + n{t) 

(3) 

with being a realization of the channel spreading 
function. In practice is causal and has finite 

support. The second order statistics of S(/i) is 



E{S(/.)SH} = C(Ai)^(/i - 



(4) 



where C(/i) is the scattering function with ||C||i = 1 
(no overall path-loss). To obtain the data symbol Xm 
the receiver does the projection 



{SAmg, r) 



iSArag)it)r{t)dt (5) 



onto a time-frequency shifted version of the function 
g{t) (i.e. an equalization filter) with ||g||2 = 1. Let 



H„ 



def 



SAmg, T-LSaui) 



(6) 



be the elements of the channel matrix H E C^^^ and 

def 

the noise n,n ~ {SAmg, n) the transmission scheme 
can be formulated as the linear equation x = Hx + n. 

2.2 Weyl-Heisenberg Signaling on 

With the connection of shift operators to unitary repre- 
sentations of the Weyl-Heisenberg group it is straight- 
forward to pass over to finite dimensional models, 
i.e. (cyclic) shift operators on which are given 
as unitary representation of finite Heisenberg groups. 
Hence, let A^(L,C) be the algebra of L x L matrices 
over C. Then the operators S^ E Ai{L,C) are given 
as the matrices 

(5p)™„- W^.e'^*^^'™"'^ (7) 

where all index-arithmetics are modulo L. The finite 
Heisenberg group is then 

Hl = k^^e^i} (8) 

where Fl = {0 . . . i - 1}. 

3 Formulation of the Problem 

In the view of multicarrier transmission only single car- 
rier equalization is considered. Interference cancellation 
is not used in this field due to complexity reasons. 
Hence it is then naturally to require a (the channel 
gain of the lattice point m E T) to be maximal and 
the interference power b from all other lattice points to 
be minimal as possible, where 



a = \Hr, 



P and &'='5] |i/™,„P (9) 

n^m 

This addresses the concept of pulse shaping, hence to 
find jointly good pulses {g, 7} (or precoders and equal- 
izers) achieving maximum channel gain and minimum 
interference power A comprehensive framework for 
the optimization of redundant precoders and equalizers 



with respect to instantaneous time-invariant channel 
realizations (assumed to be known at the transmitter) 
is given in [4]. However in certain scenarios it is much 
more realistic to adapt the pulses only to the second 
order statistics, given by C{fj,) and not to a particular 
realization This is in the sense of defining for 

the considered time-frequency slot m 



E«{a} 



(t2 + E«{&} - a2 + 5 I]-H{a} 

(10) 

as a long term performance measure. Optimal signaling 
via (1) and (3) maximizing (10) independent of S(-) 
is of central relevance for band efficient and low- 
complexity multicarrier implementation. For example, 
results on joint multipath and Doppler diversity crit- 
ically rely on so called approximate S -independent 
basis expansions proposed in [5] of Ti.. But a gen- 
eral approach how to obtain the "best basis" for Ti 
being WSSUS operators, i.e. maximizing (10), is still 
unknown. Nevertheless, some iterative methods are 
contained in [6]. 

The derivation of the lower bound in (10) can be 
found in [7] and in the context of pulse shaping in 
[2]. Equality is achieved if the set Q{'y,A,¥'^) =^ 
{SAnj}ne¥^ (called a Gabor set or family) establishes 
a tight (Gabor or Weyl-Heisenberg) frame [8], [2]. The 
constant is called the Bessel bound of (7(7, A, F^) 
and related to its redundancy. For Gi'J, A, F^) being an 
ONE it follows = 1. 

Straightforward computation yields the channel fi- 
delity (or averaged gain term) given as 



Ew{a}= / C{^,)\{g 



(11) 



and the averaged interference power 

Ew{fc} = E / C{ii)\{g,SAm+^.l)?diJL (12) 

Hence, the SINR(5,7, A) is independent of m. Let 
denote the hermitian adjoint operator of S^ with re- 
spect to (•, •). Then the channel fidelity can be rewritten 

as 



j ci^l)s,Ts;d^l 



g) "^^ TrA(r)G 



E-H{a} = (5, 

(13) 

where G (and F) is the (rank-one) orthogonal projector 

def 

onto g (and 7), i.e. Gf ~ {g, f)g. Similarly 



E«{5} = Tr 



m=iO 



Am 



Air)s 



Am 



G = TrC(r)G 



(14) 



where A{-) and C(-) are affine maps acting on linear 
operators. The definition of A{-) in particular is also 
known as Kraus representation of a completely positive 
map (see for example [9]) which establishes a relation 
to quantum channels. Due to ||C||i = 1 and 



being unitary operators, the following properties can 
be verified 

^is unital <^ A{I) =1 
^is trace preserving 4=> TrA{X) = TrX 
Ais hermiticity preserving <^ A{X*) ~ A{X)* 
Ais entropy increasing <^ A{X) ~< X 

(15) 

where -< is in the finite case the partial order due to 
eigenvalue majorization Thus, A{-) flatten the eigen- 
value distribution of its input (increasing its entropy). 
After application of A{-) onto a rank-one projector F, 
the "output in the averaged sense" (over an ensemble of 
WSSUS channels) A{r) ^ T is not (in general) rank- 
one. In this picture so called additional eigen modes 
occur which can not be collected together using a rank- 
one equalizer (a single equalization filter). 

With D{r) "^^ C(r + aH) the SINR optimization 
problem reads 

TrA(T)G 

max SINR(G,r)= max ^ (16) 
GXez.A ^ ' ' GSez,A TrD{T)G 

The maximization is performed over possible lattices A 
and G, r G Z, where Z denotes the set of orthogonal 
rank-one projectors, i.e. 

Ml = {z\Trz = l,z* = z,z>0} 

Z = {z\ze Mi,z^ = z} 

Note that the convex hull of Z is the subset Mi of 
positive-semidefinite trace class operators. 

The maximizing G for fixed F in (16) is achieved 
by an orthogonal projection onto the generahzed 
eigenspace corresponding to the maximal generalized 
eigenvalue Amax(^(r), -D(r)). Thus it remains the 
"transmitter-side only" optimization: 

max^ SINR(r) = max^ Xn...{A{r), D{r)) (17) 

With the definition of an adjoint channel it also possible 
to obtain a "receiver-side only" optimization [2]. Due 
to joint quasi-convexity of the function Ainax(-, ■) (see 
for example [10]) the constraint T G Z can be relaxed 
to the convex set Mi. Thus, the optimization problem 
is identified as convex constrained quasi-convex max- 
imization. Furthermore, if the inverse of D{r) exists, 
the problem can be rewritten as a classical eigenvalue 
problem. Note that this is non-convex optimization. 

Now, the lower bound in (10) suggests the maximiza- 
tion of the channel fidelity Tr^(r)G' only, i.e. 

max TrA(r)G = max XraUM^)) < 1 (18) 

which does not depend on the lattice A. The derivation 
from the left to the right side in (18) is again due to 
convexity of Aniax(')' linearity of A{-) and unitarity of 
S^. It can be shown that this formulation is now equiva- 
lent to the problem of maximizing the quantum channel 
fidelity [11] for G being a pure state (rank-one). Where 



the solution of (18) for single-dispersive channels is 
straightforward, the general case of this optimization 
problem - convex constrained convex maximization - is 
unsolved in general. For C((u) being a two-dimensional 
Gaussian the solution was found in [12]. Another proof 
which additional gives the uniqueness of the solution 
is in [13]. 

Already in [2] it is conjectured that with a proper 
selection of a basis for the dimensional real vector 
space of hermitian operators on the left side of (18) 
could be rewritten as a bilinear program over so called 
— 1 dimensional Bloch manifolds B^{L) [14], i.e. 



max TrA(X)Y 

X,Y&Z 



max (x, av) 



(19) 



where (a^ ) G A4{L^, M) is then the matrix representa- 
tion of A{-) in this basis. Indeed - this parameterization 
will be used in the next section. 

Finafly, if the Gabor set g{T, A, F^) with the "chan- 
nel fidelity"-optimal F would establish a tight frame 
for some lattice A, the solution maximizes SINR(F) 
too. In this case "channel fidelity"-maximization equals 
minimization of the averaged interference. Normally 
this is not the case and it remains lattice optimization 
with dct A = const. In pulse shaping procedures then 
a so called orthogonalization with respect to A has to 
be applied on F to minimize the Bessel bound B-y > 1 
[15], [2]. But it will turn out that this step is not needed 
here for the case (at det A = 1). 

4 The 2 X 2 WSSUS Channel 

For this simple toy model as the underlying 
Hilbert space is assumed, i.e. L ~ 2. The 
corresponding finite Heisenberg group is 

H2 = {5(0,0), S'(i,o),S'(o,i), 5(1,1)}, where 



'(0,0) 



^(1 



(1,0) 



1 

1 

1 

1 



'(0,1) 



(1,1) 



1 

-1 

1 

-1 



(20) 



As already intended in the introduction these matrices 
represent four basic channel operations which occur in 
a randomly weighted superposition. The matrix 5(1 0) 
represents the only cyclic shift that exists, hence it 
switches the input samples. The element 5(o,i) does the 
same in frequency domain. With 5(i.i) both operations 
occur simultaneously and 5(0,0) does not change the 
input at all. The following relations are important 



(21) 



5(0,0) — cTo 5(0,1) — 0-3 — FS(^i^o)F* 
^{1,0) = c^i -^(i.i) = *'^2 = 5(0,1)5(1,0) 

where the ai are the well known Pauli-matrices and F 
is the 2x2 Fourier matrix given as 



def 1 
^-7=2 



1 



(22) 



Furthermore we define the following constants 



Po = C(0,0) 

P2 = C(1,1) 



Pi = C(1,0) 
P3 = C(0,1) 



(23) 



given by the four possible values of scattering function 
C(/i). This allows us to write 



is a 2-sphere in R"' with origin (1, 0, 0, 0). To visualize 
why it is not easy to extent this concept to higher 
dimensions we give a short proof. 

Proof: We have to proof TrX = TrX^ = 1 and 
X > 0. 

The trace normalization represents the I2 normaliza- 
tion of the precoders and equalizers. The condition 



A{X) = J S^XS;C{^l)d^l = J2 P^<^^X< (24) TrX = i TV^ x,a, = i ^ Tra, = ^xq • 2 



The Pauli-matrices establish an orthogonal basis for the 
real vector space of hermitian 2x2 matrices with inner 
product {X, Y) = TrX*Y . Thus, every 2x2 hermitian 
matrix X has a decomposition X = ^ Si=o ^i'^i 
Xi £ M. Furthermore they establish up to factors the 
finite Weyl-Heisenberg group itself as shown in (21). 
This additional property will admit the direct solution 
of the problem. The following properties are useful to 
verify the calculations later on: 



dct (Ji 



25« 

-1 for i = 1,2,2, 
Oi j = (25) 

aj i = 



where eijk is the Levi-Civita symbol. Note that the 
Pauli-matrices are unitary and hermitian. 

4.1 The Channel FideUty 

Recall the equivalent problem formulation 



max TrA(X)Y 

x,Y&z 



^max^^(x,ay) (26) 



where B\2) = {x\^J2^ x^<J^ e Z} is a 3-dimensional 
sub-manifold in — the Bloch manifold for 2x2. The 
matrix (a,;j) S A^(4,R) is the corresponding matrix 
representation of A{-) with elements 



(27) 



The Bloch parameterization is a well known tool in 
quantum physics which admits for L ~ 2 the simple 
interpretation of 3-dimensional Bloch vectors. Because 
it is not very common in this context a short overview 
will be given. First let us evaluate the conditions for 
a vector x G R'' to be in B^{2). We will adopt the 
notation ^ which means x = {xi,X2,xs) for x — 

{xo,Xi,X2,X3) e R"*. 

Lemma 1: A real vector a; g R^ is in S^(2) iff = 
1 and ||a;||2= 1. 

Thus, it is the fact that the Bloch manifold 

B^{2) = {x = (xo, x)\xo = 1 and ||.f||2= 1} 



(28) 

is fulfilled if a;o = 1. The second trace requirement is 
the rank-one constraint if X > 0, i.e. the condition 



TrX^ = -Y,Tr^x,ai){xjn,) = \\\x\\l 



is fulfilled if ||a;||2 = 2. Using ||a;||2 — Xq 
1 + = 2 gives the requirement: 

\ml = i 



(29) 



(30) 



To ensure X > we need conditions on the determi- 
nants. Firstly 



det X = — dct Xidi 



-dct( -^0 + ^3 Xi-iX2 ^ (31) 

4 \ a;i + 1X2 - X3 



> 



which is automatically fulfilled due to xq = 1 and 
\\x\\ = 1. Secondly the upper left sub-determinant is 

i det(a;o + 2:3) = ^i^o + X3) > ^ X3 > -1 (32) 

Its non-negativity is also automatically fulfilled due to 
11x11 = 1. ■ 

Remark: The condition X > is only in the L = 2 
case automatically fulfilled. 

Next we will explicitely compute the matrix representa- 
tion (fly ) for the completely positive map A{-) in terms 
of the Pauli basis. 

Lemma 2: The matrix representation of the com- 
pletely positive map A{-) is the diagonal matrix 



.1 



1 



1 



(a^j) = diag(-,(po - -^^ (po +P2) 

(Pa+Ps) - ^) 

Proof: The matrix elements A{-) with respect to 
the Pauli basis are given as 

aki = i Tr A{<Tk)ai 



kl 



(33) 



Using the properties of the Pauli matrices one can 
compute 

Trunki = 2 < (34) 
1^+4; else 

which gives then 

{a^j) = idiag(l,po +J3i - P2 - Ps, 

Po -P1+P2 -P3, 
Po -P1-P2 +P3) 



(35) 



Using now the normaHzation, i.e. po = I —pi —P2—P3 
gives the desired result. ■ 

Theorem 1: The solution of the problem in (26) is 

max {x,ay) = \{\+ max {|2(po +]5fe) - 1|}) 

x.yeB^(2) 2 fc=l,2,3 

(36) 

Proof: Using Lemma 1 and 2 gives explicitely: 

1 > _F =^ max (x, ay) 

x,y€8^{2y 

= 7:{xoyo+ max {x,by)) 

^ l|S||2 = ||S'||2 = l 

= i(l + max{|2(po - 1|, |2(po +P2) - 1|, 

^ ^ V ' V ' 

a;(opO(i) = (i_ 1^0,0) a;<"P"(2) = (l, 0,1,0) 

|2(po+P3)-l|})>^ 
V ' 2 

2:i»P')(3) = (l, 0,0,1) 

(37) 

where (6,j) is the lower right 3x3 sub-matrix of 
a. Because (6,;j) is a diagonal matrix, only the three 
optimal vectors a;^°P'*(7i) given in (37) are possible. ■ 

Furthermore the fcth component of the optimal equal- 
izer is 

y'°P''(n)fe ^ sign{2(po + Pn) - \}x^°^'\n)u (38) 
4.1.1 Discussion 

The simple solution (37) is well-suited to discuss 
several cases which occur also in the general WSSUS 
pulse shaping problem. 

(1) non-dispersive channels: If p^ ~ 1 for some 
k = ... 3, this will yield F = 1 which is achieved 
with any This is the so called flat fading 
(non-selective) channel. The maximal channel fidelity 
("F = 1") can be achieved. No precoding is needed. 

(2) single-dispersive channels: If pk = pi = for some 
k,l = ... 3 and k 7^ / yields again F = 1 which 
achieved with a;*°P''(ri) where n ^ k and n ^ I. This 
is for example for k ~ 2 and I = 3 the frequency 
selective, time-invariant channel. All contributions S'^ 
can be diagonalized simultaneously which is achieved 
for example in OFDM. The maximal channel fidelity 
("F = 1") is achieved again (in practice there is still a 



fidelity loss for OFDM due to the cyclic prefix). 

(3) doubly-dispersive "underspread" channels: If in 
general only Pk > ^ for some k ~ ... 3 we have 
only F > i. But if fc = follows F = po + pn where 
n = argmax,„=i,2,3{pm}, achieved with a;^°P'^(n). A 
closer inspection shows, that if fc 7^ the solution will 
depend on some ordering property of the scattering 
powers. The optimal precoder depends now explicity 
on multiple values p,„. 

(4) Completely overspread channels: If pk ^ \ for all 
k the channel fidelity gives F = i which is achieved 
with any This is the worst case scenario. 

It is quite interesting what happens if we fix the 
scattering power po> i-e. to consider F as the function 
F{po,p) where p = {pi,P2,P3) 

F{po,p) = 1(1 + max {\2{po - Pk) ~ 1|}) (39) 

Clearly, F{po,p) is jointly and separately convex in po 
and p. Furthermore is F{po,p) = F{po, Hp) for every 
permutation 11, so F is Schur-convex in the second 
argument (see for example [16]), i.e. for a fixed po 
follows 



Pl>- P2^ F{po,pi) > F(po,P2) 



(40) 



Using Schur-convexity for every fixed po follows: 
(5) the worst case channel: is given for Pk = for 
all fc = 1, 2, 3 yielding 



mill F(po,p) = i^(po,^^(l,l,l)) 

|p||i = l-po o 

1. 1, ,N 1 2, 1, 



(41) 



Note that in the quantum context this corresponds to the 
general depolarizing channel or Lie algebra channel 
[14]. 

(6) the best case channel: is given for pk ~ (l — po)ek 
where Ck is a standard basis vector and k E {1,2,3} 
yielding 

max F{po,p) = F{po, (1 - po)4) 
IIPIU=1-P0 ^ ^^2) 

= -(l + max{|2po-l|,l}) = l 

achieved with .T'°P''(fc). This is again an single- 
dispersive channel because <tq commutes with ak- 

4.1.2 Construction of the precoders 

In this part we will explicitely calculate the correspond- 
ing matrix representation X'°P''(n) from the Bloch pa- 
rameterizations a;*°P''(ri) using X = ^^i^i'^i which 
gives 



X(°P')(„) = -(ao + a„) 



(43) 



These matrices are the rank-one projectors onto the op- 
timal precoders in the original problem. Hence, turning 
back to precoders and equalizers x{n),y{n) G it 
can be verified that similarly) 

*"-'"">^;^(");^(")* 

^(opt).oN 1 ( 1/2(1 Tl) \ 1 ( 1/2(1 Tl) V 
^ V2(l±l) 1/2(1±1) ) 

(44) 

In the following I will use w.l.o.g. the i-version, hence 
X*°P''(n) = y'^°P''(n). The following precoders are then 
solutions of the optimization problem 




The solution a;(l) is maximally localized in the fre- 
quency domain, i.e. completely spread out in the time 
domain. The reverse case holds for .t(3). To understand 
x{2) one has to perform a rotation in the time-frequency 
plane. 

4.2 The Multiplexing Scheme 

Solving the channel fidelity problem as done in the 
previous section is the first step toward optimal sig- 
natures for a given WSSUS channel statistics. In this 
part we will discuss now how multiplexing has be per- 
formed. In C^-case at spectral efficiency dctA^^ = 1 
this means multiplexing of a second data stream only. 
Recall that with Gabor (Weyl-Heisenberg) signaling 
the same channel fidelity is achieved for the second 
data stream. It remains first to select the lattice with 
minimum averaged interference and then perform some 
kind of orthogonalization ("tighten") procedure. But the 
case admits already the following transmitter side 
orthogonality relations 

(a;(n), S'(o,i)x(n)) = i Trcr3(cro + (^n) = <5n3 

{x{n),S(i^o)x{n)) = i Trcri(cro + o-ri) = <5>ii (46) 

{x(n),S(x,i)x{n)) = - Trcr2((To + o-„) = iSn2 

Hence, for each channel optimal pulse "n" there are 
several schemes "n," for multiplexing a second data 
stream which admits transmitter-side orthogonality, i.e. 
no further orthogonalization is required. The "channel 
fidelity"-optimal precoder is also SINR-optimal. If for 
example the channel optimal precoder is x(3), "time- 
division multiplexing" via 5'(o.i) is one of the optimal 
schemes. For a;(l) in turn "frequency-division multi- 
plexing" is the right scheme. 



5 Conclusions 

In this article new insight into the WSSUS pulse 
shaping problem are given from a precoding viewpoint. 
The precoding problem is solved for a very simple 
class of random 2x2 channels under the assumption 
that the transmitter has only knowledge of the second 
order statistics and the receiver has full knowledge 
of the channel. It is observed that optimality in the 
channel-fidelity sense is achieved with concentration in 
a certain domain similarly as one would expect from the 
continuous case. Unfortunately the direct extension of 
this approach to higher dimension is problematic due to 
the difficulties arising with the Bloch parameterization. 
But a staggered extension could be conceivably and will 
probably studied in future. 
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